Relevance of many-body interactions for correlated electrons 
in the strong-coupling limit 
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Many-body interactions in effective field theories for disordered interacting electrons are consid- 
ered. It is shown that three-body and higher interaction terms are generated in perturbation theory, 
and some of the physical consequences of these interactions are discussed. It is shown in particu- 
lar that they will in general be important for any effects governed by strong-coupling fixed points. 
This implies that the usual generalized nonlinear sigma-model for disordered electron systems is 
incomplete, and not suitable for studying strong-coupling effects. 
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INTRODUCTION 



Effective field theories are a very useful tool, both in 
statistical mechanics and in particle physicsi The basic 
idea is to construct an "effective" theory, valid at large 
length scales and long times, which contains only those 
degrees that are important in this regime, while all oth- 
ers have been integrated out. Since the effective theory 
is simpler than the underlying fundamental theory that 
contains all degrees of freedom explicitly, it is easier to 
solve. If the fundamental theory is known, as is the case 
in condensed matter physics, effective theories can be de- 
rived from it. If it is not known, as in particle physics, 
the effective theory can be guessed with feedback from 
experiments. In either case the effective theory in general 
contains features that are not present in the fundamen- 
tal or microscopic theory. For instance, interaction events 
between particles that are sequential occurences of funda- 
mental interactions on microscopic scales will appear as 
basic interactions on the coarse-grained length and time 
scales of the effective theory. An example from particle 
physics is the Fermi theory of beta decay, which assumed 
a point-like interaction between the particles involved. 2 
Later, in more microscopic theories of the weak inter- 
action, it became clear that there is internal structure 
in Fermi's interaction related to the exchange of gauge 
bosons £ 

In condensed matter physics, the only interaction in 
the microscopic theory is the Coulomb interaction. We 
will be concerned with electron-electron interactions in 
disordered metals, and therefore we take the "fundamen- 
tal" interaction to be the screened Coulomb interaction. 
Let us consider processes in which two electrons interact 
at some point in space and time, and some time later 
a third electron interacts with one of the two some dis- 
tance away from the first interaction point. In an effec- 
tive theory that has integrated out the behavior at short 
length and time scales, such a process will appear as a 



"fundamental" interaction between three electrons, since 
the effective theory can no longer resolve the individual 
microscopic interaction processes. In classical statistical 
mechanics the importance of such effective many-body 
interactions is well known. An example is the expansion 
of transport coefficients in powers of the particle number 
density. To obtain the contribution at any given (suffi- 
ciently high) order in the density one needs to consider 
collisions between arbitrarily many particles^ Analogous 
effects have been considered for many-electron systems^ 
although the connection with effective many-body inter- 
actions was not made explicit. Furthermore, the con- 
struction of a complete effective theory requires that any 
many-body interactions that are generated in perturba- 
tion theory be included in further iterations of the renor- 
malization process that integrates out the short-range 
degrees of freedom. This has never been done; existing 
effective theories for disordered many-electron systems 
contain two-body interactions onlypl 

In the present paper, we show explicitly that many- 
body interactions in such systems are generated under 
renormalization. The many-body interactions generated 
are of long range in space and time due to the diffusive 
electron dynamics. As a consequence of their long-range 
nature, the naive renormalization-group (RG) scale di- 
mensions of these terms vanish in two-dimensions, which 
implies that they should be important in theories of 
the Anderson-Mott metal-insulator transition near two- 
dimensions & We will clarify in what sense this is the case. 
We will further show that these many-body interactions 
can lead to qualitatively new scaling behavior if the in- 
teractions are strong enough. 

The organization of this paper is as follows. In Sec. 
ITT1 we give simple physical arguments for the existence 
of effective many-body interactions and their expected 
structure. In Sec. IIIII we perform an explicit calculation 
showing that such terms are indeed generated in pertur- 
bation theory, starting with a model that has two-body 
interactions only. In Sec. lIVI we discuss the physical con- 
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sequences of these terms, and in particular their relevance 
for strong-coupling problems. Section IVI contains a con- 
clusion, and in the appendix we discuss some aspects of 
</> 4 -theory that are analogous in some respects to our per- 
turbative calculation. 



II. PHYSICAL ARGUMENTS FOR THE 
EXISTENCE OF EFFECTIVE MANY-BODY 
INTERACTIONS 

Let us consider an ensemble of interacting electrons 
in the presence of quenched disorder. For simplicity, we 
will model the screened Coulomb interaction by an in- 
stantaneous, point-like model interaction whose coupling 
constant we denote by K^ 2 \ The action will therefore 
contain a term 



<j(2) ^ 



J dxdy J dr n{x,r) S(x — y)n{y 1 r) 
K {2) J dxdy Tj2 n ( x > tt n ) 6{x - y) n{y, -Q r , 



(2.1) 



Here n is the electron number density field, which is 
a function of position x and imaginary time r, and 
f3 = l/T is the inverse temperature. We use units 
such that ti = kj} = 1. In the second line we have 
performed a Fourier transform from imaginary time to 
bosonic Matubara frequencies f2 n = 2itTn. 

If neither the disorder nor the interaction is too 
strong, 19 the dynamics of the electrons will be diffusive. 
This means there are particle-hole excitations, or diffu- 
sons, that are described by a diffusive propagator 



V n (x -y) = 5{x - y) {-DV 2 + \Q n \y 



(2.2) 



The exchange of diffusons then provides an effective long- 
ranged interaction between the electrons. Consider, for 
instance, three electrons that are pairwise coupled by dif- 
fusion propagators, see Fig.Q] Each two-body interaction 
carries an amplitude K^ 2 \ and we therefore expect this 
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FIG. 1: An effective three-body interaction mediated by three 
diffusion propagators. 
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FIG. 2: An effective three-body interaction mediated by two 
diffusion propagators. 



process to contribute a term to the effective action, 



(34) 



h )' J dx dy dzj^drTY^ 



n(x, t) 



xV m (x - y) n(y, r) V m (y - z) n(z, r) 
xV rn (z-x) . (2.3) 

Here we have localized the diffusons in imaginary time, 
i.e., we have neglected their dependence on the external 
frequency arguments. This implies an effective three- 
electron interaction amplitude, defined in analogy to Eq. 
(12.1(1 . that is given by 

K^{x,y,z) = (K^yT^2v m (x-y)V m (y-z) 

m 

xV m (z-x) (2.4a) 

In an effective theory that cannot resolve the positions x, 
y, and z, i^ 3 ' 1 ) will appear as a point-like three-electron 
interaction. For later referece we note that in momentum 
space, and with the external momenta put equal to zero, 
i^ 3,1 ) reads 

^) = (irW) 3 T^I^(P m (k)) 3 . (2.4b) 

m k 

Notice that the frequency-momentum integral in this ex- 
pression is infrared divergent in all spatial dimensions 
d < 4. This singularity will be cut off by any nonzero ex- 
ternal momenta and frequencies. Physically, this means 
that the three-body interaction is of long range in space 
and time, as was mentioned in the Introduction. We will 
come back to this point in Sees. Illll and IIVI below. 

While i^ 3 ' 1 ) is perhaps the most obvious three-body 
interaction term, it is easy to see that there are others, 
including terms that are only quadratic in the two-body 
interaction amplitude K^ 2 ' . Consider the situation in 
Fig. [21 where two electrons coupled by the original short- 
ranged two-body interaction interact with a third one 
by exchanging diffusons. Since there are two electrons 
at the same point in space in this process, we expect 
the long-ranged effective interaction to be mediated by a 
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diffuson squared. This leads to an effective three-electron 
interaction amplitude, at zero external frequencies, 

K^\x,y,z)= (k^) 2 5{x-y) (V n=Q (y - z)f . 

(2.5) 

It is obvious that there cannot be any three-electron in- 
teraction terms that are linear in K^ 2 \ 

Analogous arguments lead to the conclusion that there 
are four-electron interaction terms, starting at order 

(ii^ 2 )) , etc. In the following section we will ascertain 
the existence of such many-body interactions by means of 
an explicit perturbative calculation for a specific model. 



III. GENERATION OF MANY-BODY 
INTERACTIONS IN PERTURBATION THEORY 

A. Generalized nonlinear sigma-model 

We now turn to an explicit calculation that shows 
how many-body interactions are generated by renormal- 
izing models that contain two-body interactions only. 
For definiteness, we take as our starting point the gen- 
eralized nonlinear sigma-model for disordered, interact- 
ing electrons! that has been used extensively to de- 
scribe metal-insulator transitions as well as magnetic 
transitions 6 in solids. The action reads, 



dx tr 



(vQ(x] 



2H / dx tr 



(3.1a) 



Here Q is a hermitian matrix field subject to the con- 
straints 



Q 2 {x) = 1 , trQ(x)=0 . (3.1b) 

The matrix elements Q"m carry fermionic Matsubara fre- 
quency labels n, m, and replica labels a to deal with the 
quenched disorder. The matrix elements arc themselves 
four-by-four matrices to allow for spin and particle-hole 
degrees of freedom. They constitute the soft-mode com- 
ponents of an underlying matrix field Q that comprises 
bilinear products of fermionic fields ip and tp according 
to the correspondence 

/-V'itV^t -V'ltV'aj. V'itV'2t\ 

_ i -V'uV^t -iptifoi -ipnfai ^uV'at 
( " 12 ~~ 2 ^l|-02T -011^21 dufol --011-02! 

. --0l1-02t --01T-021 ~1pl-[1p2l 01T02T 

(3.2) 

Here all fields are understood to be taken at position 
x, and 1 = (m,a!i), etc., comprises both frequency and 
replica labels. It is convenient to expand the 4x4 ma- 
trices in a spin-quaternion basis, 



/ 



o 

?12(x) = ^ ( T r 



SV2 



(x) 



(3.3) 



Here To = sq = I2 is the 2x2 unit matrix, and 
Tj = —Sj = —iuj, (j = 1,2,3), with 01,2,3 the Pauli 
matrices. In this basis, i = and i = 1,2,3 describe 
the spin-singlet and the spin-triplet, respectively. An ex- 
plicit calculation reveals that r = 0, 3 corresponds to the 
particle-hole channel (i.e., products tpip), while r = 1,2 
describes the particle-particle channel (i.e., products iptp 
or ipifi). In this basis, the electron number density field 
as a function of x and a bosonic Matsubara frequency 
Q n is represented by 



r=0,3 



{V^iy Vtr (r r (g)s ) Q 



m,m+n \ 



(3.4) 

G and H in Eq. H3.1afl are coupling constants that rep- 
resent the disorder strength and the frequency coupling, 
respectively. Their bare values are proportional to the re- 
sistivity in Boltzmann approximation, and to the density 
of states in self-consistent Born approximation, respec- 
tively. O is a frequency matrix with matrix elements 



^12 = (to ® s ) S 12 2nT{ ni + 1/2) 



(3.5) 



The final term in Eq. (|3.1a|l describes the two-body 
electron-electron interaction. From Eqs. (|3.2I3.4|I it is 

clear that it must be quadratic in Q. If one separates the 
interaction into a spin-singlet interaction between num- 
ber densities, and a spin-triplet interaction between spin 
densities, A- ln t reads 



with 



A 



(2,s) 



A 



ttT 



(2) _ - (2,s) 



A ( ' ' 4- A 

Ant I 



(2,t) 



(3.6a) 



#(2,s) 



; E 

r=0,3 



E E 



ni,Ti2,m a 



A 



(2,t) 



tr (T r ®«o)QST, ni+w W 



tr ((T r ^o)g: +min2 (x))] , (3.6b) 



f ^)/ rfx ^(-ir E EE 

r— 0,3 Tii,ri2,m a i—1 

((rrOsi)g^„ 1+m (x)' 



tr 



tr ((t p <8« 4 )Q™ + „ a (x) 



(3.6c) 



Here if( 2 > s ) and _ftT( 2 > s ) are the spin-singlet and spin- 
triplet two-body interaction amplitudes, respectively. In 
general, there also is an interaction amplitude in the 
particle-particle channel, which we neglect here. 

Finally, for explicit calculations it is convenient to elim- 
inate the constraints given by Eq. Ij3.1b|) . This can be 
done by means of the block matrix parametrization 



r,i— 



Q 



(3.7) 
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Here the four block matrices represent, clockwise from 
upper left, the matrix elements of Q with frequency labels 
n\,n-2 > 0, n\ > 0, n-i < 0, ni, n 2 < 0, and ri\ < 0, n 2 > 
0. 

From Eq. (|3.4|l it follows that a point-like, instanta- 
neous three-body interaction term involving three num- 
ber density fluctuations would take the form 

43 = ^^ s) h* E (^)" +s+t E 

J r,s.t=Q.3 "i,"2."3 

x<5„ 1+ „ 3+ „ 5 ,„ 2+ „ 4+ „ 6 [ tr ( T r ® so) Qm„ 2 (aO] 

a 

x [tr (t s <g> s ) Q««„ 4 (x)] [tr (r t ® s ) Q^„ 6 0*0] . 

(3.8) 

We will now show that such a term is indeed produced 
by renormalizing the bare action A given in Eq. (|3.1a(l . 




FIG. 3: Diagrams that renormalize the two-body interaction. 




FIG. 4: Diagrams that can generate an effective three-body 
interaction. 



B. Loop expansion 

To proceed, we expand the action A, Eq. I|3.1a[l . in 
powers of q. To Gaussian order we obtain a quadratic 
form whose inverse determines the Gaussian propagators. 
In Fourier space, the latter read 



(r9i 2 (Pl)iQ , 34(P2)) = -j 
with 



G i 
-S rs 8 l3 l r M~l M ( Pl ) 



(3.9a) 



5i 3 V ni - n2 (p) 



x V: 



(p) 



-S13S24, I'm -712 (p) 



- n 2 
(3.9b) 

(3.9c) 



Here v = s, 1/1,2,3 
propagators 



t, and we have introduced the 



V n (p) = l/(p 2 +GHfl n ) , 
V s n (p) = l/(p 2 +G(H + K^)n, 

2>*(p) = l/(p 2 +G(H + K^)n r 



(3.9d) 
(3.9e) 

(3.9f) 



which are proportional to the basic diffusion propagator 
defined in Eq. 

We now perform a systematic loop expansion, and con- 
centrate on the renormalizations of the interaction terms. 
For simplicity, we consider only the particle-hole chan- 
nel degrees of freedom, i.e., we neglect the propagators 
with r = 1,2 above. A physical situation that realizes 
this approximation is, e.g., a system with magnetic im- 
purities, which give the particle-particle channel propa- 
gators a mass, so they drop out of the soft-mode effective 
theoryi^ 



At one-loop order, the two-body interactions get renor- 
malized, and they acquire frequency and momentum de- 
pendences in the process. This effect is due to the dia- 
grams shown in Fig. |3J and it has been studied in detail 
before^ However, there also are three-body interactions 
generated in the process. Let us concentrate on the pure 
spin-singlet term given by Eq. I|3.8|) . To lowest order in 
powers of q, it will manifest itself in particular in a term 

(AttT) 2 22 $nx+n 3 -n 5 ,n 7 .+n A -n fi J dx dy dz 

114,715,715 

X ^n 1 'nln 3 ,n 4 n 5 n 6 { X l Vl Z ) a1 ni n 2 ( X ) oQn 3 n 4 ClO 0#n 5 n 6 ( Z ) L 

(3.10) 

Here we have allowed for a frequency and real-space de- 
pendence of the three-body interaction amplitude, and 
for simplicity we only consider the ®q components of the 
matrices q. A vertex with the structure of Eq. I|3.10|l can 
in principle be generated by any of the three diagrams 
shown in Fig. 0] 

It is easy to see that diagram (a) in Fig. 0] does not 
contribute to a three-body interaction; the realizations 
of this diagram that have the correct replica structure 
do not have a frequency structure consistent with Eq. 
(|3.10l) . This is consistent with our conclusion, in Sec. ITTl 
that there are no contributions to that are linear 

in K^ 2 \ Of the remaining two diagrams, (c) is at least 
of cubic order in K^ 2 \ and in particular produces terms 
that have the structure of i^ 3 ' 1 ) in Eq. I|2.4a|l . Diagram 
(b) has contributions of the correct structure that are 
of 0((K( 2 >) 2 ), as well as contributions of higher order. 
The easiest check for the existence of therefore con- 
sists of a systematic calculation of diagram (b), keeping 
only terms of second order in the two-body interaction 
amplitude K^ 2 ' s \ Such a calculation yields 
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^123,456 ( X ' V' Z ) = K U3A56( X >yi Z ) + K Ull56(y> X > Z ) (3.11a) 

where K 1 - 3 '^ reads, in Fourier space, 

#123,456(kl,k 2 ,k 3 ) = (GK s /8) 2 f 

Jp 

x {-26(5 > 1 - 4) [(p + kO 2 + (p + k 2 ) 2 + G#Qi_ 2 + GH^ 2 ] (2? 3 - 2 (p)) 2 £>i- 4 (p - k 3 ) 

+26(5 < 1 - 4) [(p + kO 2 + (p + k 2 ) 2 + GHQx.i + dm 3 _ 4 ] (2?i- 4 (p)) 2 £> 3 - 2 (p - k 3 ) 
+26(5 > 1 + 3) [2? B -i(p) 2?i_ 6 (p - k 3 ) + £> 5 - 3 (p) £>5- 4 (p - k 2 )] 
-26(5 < 1 + 3) [X>b_ 2 (p) X>2- 6 (P - k 3 ) + Ps-6(p) ^4-e(p - k 2 )] 
-26(5 > 3 - 2) [X>6_ 3 (p) £>3-e(p - k 3 ) + P 5 - 3 (p) £>s-4(p - k 2 )] 
+26(5 < 3 - 2) [X>b_ 4 (p) X>4- 6 (P - k 3 ) + £>3-e(p) 2?4-e(p - k 2 )] 

-26(5 > 1 - 2) 2?8-4(p) ^5-4-i+ 2 (p - k t ) + 26(5 < 1 - 2) X>i_ fl (p) Pi- 6 - 3+4 (p - k 2 ) 
+2P 5 - 4 (p)I?5-4+i- 2 (p + ki) - 2Pi_ 6 (p)Pi_ 6+3 _ 4 (p + k 2 )} . (3.11b) 



Here J = J dp/(2ir) d , 1 = m, etc., and the symbols 
6(5 > 1 — 4) = <d(n 5 — ni + n±), etc., with the second 
6 denoting the usual Heavyside step function, express 
constraints among the frequencies. 

This result demonstrates that a term with the struc- 
ture of K^ 3 ' 2 \ Eq. (|2.5|l . gets indeed generated upon 
renormalization of an action with a pure two-body in- 
teraction. In addition, there exist terms that repre- 
sent K { ^\ Eq. J233), as well as spin-triplet and mixed 
singlet-triplet three-body interactions. It is also plausible 
that four- and higher-body interaction terms are gener- 
ated by the same mechanism, and the existence of par- 
ticular diagrams with the appropriate structure is easily 
verified. 

We note that the momentum integral in Eq. (|3.11b|) 
diverges for small external wavevectors k, or small exter- 
nal frequencies f2, as 1/fe 2 or I/O, in agreement with the 
remark after Eq. (|2.4b|) . The RG interpretation of this 
divergence is given in the next section. 



IV. PHYSICAL EFFECTS DUE TO EFFECTIVE 
MANY-BODY INTERACTIONS 

A. Structure of renormalization-group flow 
equations 

Since structurally new terms have appeared in our ac- 
tion under renormalization, we need to add these terms 
to the action and start the renormalization process over 
again. The action as given by Eqs. 1)3. 113. 6|) thus must 
be augmented by Eq. I|3.8|l . 20 The renormalization of this 
action then proceeds along standard lines. The result is 
obviously a generalization of the known flow equations 
for the model with two-particle interactions only. For 
our present purposes we are interested only in the general 



structure of these flow equations, which can be obtained 
without a detailed calculation. 

We choose the scale dimension of a length L to be 
[L] = 1, and that of imaginary time r to be [r] = d in 
d spatial dimensions^ The field q(x) we choose to be 
dimensionless. The bare scale dimension of G is then 
d — 2 = e, the bare scale dimensions of H, K^ 2 ^ and 
are zero. The bare scale dimension of is — d, 
due to the extra factor of T that appears in the three- 
body interaction term, Eq. (|3.8|l . compared to the two- 
body interaction. If we denote the renormalized, scale 
dependent counterparts of these coupling constants by g, 
h, k s , kt, and fc 3 , respectively, we thus have to zero-loop 
order 



dg_ 

dl 
dh 

17 

dk 3 
dl 



dk s dkt 
~dl^~dl 

-(2 + e)fc 3 



(4.1a) 
(4.1b) 
(4.1c) 



Here I = In o, with b the renormalization-group length 
rescaling factor. 

To find the higher-loop order terms explicitly requires 
a detailed calculation. For k% = 0, the result is known 
completely to one-loop order, and selectively to two-loop 
order £ For all universality classes, the structure is, 



dg 

dl 
dh 

~dl 

dl 
dl 



-£5 + 5 2 /i 1) (7 s ,7t) +5 3 /f (7s,7t) ,(4.2a) 
^/ / ( l 1) (7s,7t) + % 2 // l 2) (7s,7t) , (4.2b) 
3/ s (1) (7s,7t)+S 2 /i 2) (7 s ,7t) , (4.2c) 
5/t (1) (7 s ,7t)+5 2 /t (2) (7 s ,7t) , (4.2d) 
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where 7„t = k St t/h. We note that —1 < 7 S < 0, and 
7t > 0, and the various functions / are well-behaved in 
the limit 7 S t — > 0£ 

In the presence of ka, we need to consider, first, the 
influence of k% on the flow of the other coupling con- 
stants, and, second, the fc3-flow equation itself. Simple 
counting arguments show that k$ cannot produce sin- 
gular (in d = 2) one-loop renormalizations of the other 
coupling constants. For instance, consider the second di- 
agram in Fig. with one of the vertices replaced by a 
three-body interaction. Due to the additional factor of 
T in Eq. I|3.8|) compared to Eqs. (|3.6[l this diagram will 
have an extra frequency integration compared to the dia- 
gram with both vertices given by two-body interactions, 
and will thus not be infrared singular. The structure of 
the flow equation for k% itself is therefore more impor- 
tant than the modifications of Eqs. H4.2al - IT.2dJl . and the 
crucial question is whether is it possible to overcome the 
negative bare scale dimension of fc 3 . The most interest- 
ing term is therefore the one-loop renormalization of k^ 
that is proportional to k 3 itself. Such terms exist; they 
are realized, e.g., by diagrams (b) and (c) in Fig. 0] with 
one of the three-point vertices replaced by a three-body 
interaction. Simple counting arguments show that the 
structure of the fc3-flow equation to one-loop order is 

= "(2 + e)h + k 3 gf^( % , 7t) + g 2 fP(j B , Jt,h) . 

(4.2e) 

Note that the last term on the right-hand side of Eq. 
(|4.2e|l is independent of k%. It represents the contribu- 
tions that generate £; 3 in the first place, for instance, the 
one given by Eqs. 

In general, adding a new RG variable to a set of flow 
equations can have any one of three effects by virtue of 
the new eigenvalue it adds to the set of equations lin- 
earized about any fixed point. First, it may be truly ir- 
relevant in the sense that it does not qualitatively change 
any aspects of the RG flow in its absence. Second, it may 
be irrelevant with respect to a fixed point that exists in 
its absence, but change the flow outside of the basin of at- 
traction of this fixed point. (It will in general also change 
the size of this basin of attraction.) Third, it may be rele- 
vant with respect to the original fixed point. In the latter 
two cases, outside of the basin of attraction of the origi- 
nal fixed point, if any, it may either lead to a new fixed 
point, or to flow towards strong coupling. With this in 
mind, we next discuss possible types of fixed points of 
the above flow equations. 

B. Weak-coupling fixed points 

In the usual perturbative RG treatment one looks for 
fixed points of the flow equations, Eqs. (|4.2I) . where g is 
small of 0(e), and 7 S and 7t are at most of 0(1). This 
is our definition of a weak-coupling fixed points It fol- 
lows from Eq. (|4.2e|l that the new scaling operator intro- 
duced by the presence of k^ will have a scale dimension 



of — 2 + O(e) with respect to such a fixed point. This 
is assured by the bare scale dimension of A; 3 , [/c 3 ] = — d, 
which cannot be overcome by the small one-loop term. 
In this context it is important to mention that k^ itself 
does have a component that is marginal in d — 2. This 
follows from the fact that the one-loop term in Eq. 14.2e|) 
has a contribution that is independent of k^, see the re- 
mark after that equation, and Eqs. (|3 - 1 1 1) . However, this 
component just reflects the scaling behavior of the other 
coupling constants, g, h, k s , and kt, and it does not lead 
to a new eigenvalue of the linearized RG flow equations. 
In other words, k% is not a proper scaling operator, and 
the scaling operator related to k^ has the components 
that are marginal in d — 2 projected out. An analogous 
phenomenon in 4 -theory is discussed in the Appendix. 

We also note that, alternatively, one could treat the 
three-body interaction generated by Eqs. (|3 - 1 1 1) as a truly 
long-ranged interaction with a bare scale dimension of 
— 2e. Such a procedure would lead to the same conclu- 
sion, since the part that is marginal in d — 2 would just 
reflect the scaling behavior of the two-body interaction 
amplitudes. The one-loop term independent of k% thus 
reflects the long-range nature of the RG-generated three- 
body interaction; see also the remark at the end of Sec. 
IIIII This observation justifies our using a point-like three- 
body interaction amplitude despite the fact that the one 
generated in perturbation theory was long-ranged. In the 
Appendix we discuss a similar feature of </> -theory. 

We conclude that the many-body interactions are in- 
deed important for weak-coupling fixed points, as one 
would have expected. However, since they do not lead 
to new marginal (in d = 2) scaling operators, the rele- 
vant physics is already contained in the renormalization 
of the two-body interaction constants. Weak-coupling 
fixed points will thus always be perturbatively stable with 
respect to k%, and also with respect to higher many-body 
interactions, i.e., they have a finite basin of attraction. 
This in turn implies that the coupling of k$ to the other 
coupling constants cannot change the critical behavior, 
it will merely lead to power-law corrections to scaling. In 
particular, all of the perturbative metal-insulator transi- 
tion fixed points that are known to exist for the general- 
ized nonlinear sigma- model defined in Sec. IIIII are stable 
with respect to k 3 . 



C. Strong-coupling fixed points 

Let us now consider strong-coupling fixed points, 
where the ratio of successive terms in the loop expan- 
sion is not necessarily some power of e. This can hap- 
pen if the fixed point value of j t is large, of 0(l/e), or 
infinite, even if g is still of O(e). Of course, another pos- 
sibility is that g = O(l). No controlled theories exist 
of metal-insulator transitions that are governed by such 
fixed points, but they are structurally clearly possible. 
Explicit, if uncontrolled, examples of fixed points where 
both g and an interaction coupling constant are of O(l) 
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have been given in Refs. ^3 an d It is also believed 
that strong-coupling physics governs the behavior in cer- 
tain 2-d sytems, where the interaction strength may be 
the dominant energy scale in the problem^ 

The arguments in the previous subsection that en- 
sure the irrelevance of ks obviously break down for such 
strong-coupling fixed points. We stress that this may 
be true even if the fixed-point value of g is still small. 
The point is that with, e.g., g — 0(e) and 7 t = 0(l/e), 
<77t = 0(1), and hence the one-loop term in Eq. 14.2e|) 
can overwhelm the bare scale dimension of ^3. The same 
arguments hold for the higher many-body interaction 
terms, although to a lesser degree, since for them a larger 
negative bare scale dimension must be overcome in order 
to make them relevant. 

We conclude that the many-body interaction terms 
cannot be dismissed a priori in any strong-coupling 
regime, where the dimensionless interaction amplitudes 
are large, even if the disorder is still small. In particu- 
lar, they are likely to play a role in the resolution of the 
two-dimensional metal-insulator transition problem. 



V. SUMMARY, AND CONCLUSION 

To summarize, we have shown that many-body interac- 
tions are generated under renormalization of an action for 
interacting disordered electrons that contains two-body 
interactions only. Such interactions turn out to be irrel- 
evant with respect to the perturbative fixed points that 
describe metal-insulator transitions in d = 2 + e. How- 
ever, they need to be examinated, and they likely con- 
tribute to the leading behavior, in any strong-coupling 
theory. This implies in particular that even if strong- 
coupling solutions for the generalized nonlinear sigma- 
model with two-body interactions could be found, such 
solutions would be incomplete, and probably physically 
wrong. The task of determining the strong-coupling be- 
havior of such systems, and in particular the situation in 
d = 2, is thus even harder than previously assumed. 

We conclude by recapitulating two aspects of our tech- 
nical development that can easily lead to confusion. First 
of all, let us come back to the infrared divergence of the 
integral in Eq. (|3.11b|l . Naively, this infrared divergence 
seems to offset the extra factor of temperature compared 
to the two-body interaction, making the three-body in- 
teraction marginal by power counting in d = 2. The 
same argument applies to higher many-body interaction 
amplitudes, which carry additional factors of tempera- 
ture, but come with even more divergent loop integrals. 
As we show explicitly by analyzing an analogous effect in 
</> 4 -theory in the appendix, this simple argument is falla- 
cious and the many-body interactions are perturbatively 
irrelevant, but they are likely to play an important role 
in nonperturbative regimes. Second, we have focussed 
on one particular three-body interaction term, namely, 
a spin-singlet three-body interaction. For this term, we 
have calculated all contributions to second order in the 



spin-singlet two-body interaction within a well-defined 
model. This proves the existence of many-body inter- 
actions in effective field theories for electrons, but our 
calculation is sensitive to only a small subclass of such 
terms. Many-body interactions coupling four and more 
electrons, spin-triplet interactions, and terms coupling 
singlet and triplet density fluctuations certainly exist, 
and they all need to be examined in order to system- 
atically deal with strong-coupling effects. 

Finally, we note that existing theories of the Anderson- 
Mott transition seem to lead to the conclusion that it is 
very similar to an Anderson transition, 2 - However, it is 
reasonable to assume that in strongly correlated systems 
the metal-insulator transition should more closely resem- 
ble a Mott transition^ This suggests the existence of 
a saddle in parameter space separating the fixed points 
that describe the two transitions. The disordered Mott 
transition is likely described by a strong-coupling fixed 
point. As discussed after Eqs. (|4.2|l . on the strong- 
coupling side of such a saddle the many-particle interac- 
tions are likely to be important, even if the weak-coupling 
Anderson or Anderson-Mott fixed point is locally stable 
with respect to them. In this context it is interesting to 
note that a Landau theory for the Anderson-Mott tran- 
sition in high dimensions (d > 6)MJ£ found indeed that 
for weak interactions, an Anderson transition takes place 
with increasing disorder, while for strong interactions the 
metal-insulator transition has a different nature. It is 
likely that the many-body interactions discussed in this 
paper are important for understanding the missing link 
between this theory in high dimensions, and the usual 
treatment of the metal-insulator transition problem near 
two-dimensions. 
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APPENDIX A: GENERATION OF 
HIGHER-ORDER TERMS IN 4 -THEORY 

In this appendix we recall some aspects of (/> 4 -theory 
that are analogous to the generation of many-body in- 
terac tions d iscussed in the main part of the paper. See 
Refs. llfill7llill for a derivation of the results summarized 
below. 

Consider scalar <^) 4 -theory, with an action 




(Al) 

Upon renormalization, a </> 6 -term with a coupling con- 
stant Uq is generated by means of the diagram shown 
in Fig. [SJ This diagram represents a momentum inte- 



8 




FIG. 5: Diagram that generates a </> 6 -term from 4 -terms. 



gral over three propagators. In perturbation theory, and 
at criticality, this integral is strongly infrared divergent. 
Naively, this raises the question whether u 6 can really be 
irrelevant with respect to the Wilson-Fisher fixed point. 
To investigate this, we consider the RG flow equations. 
Adding a </> 6 -term to the action, one easily finds, within 
a momentum-shell RG, and to one-loop order, 



dr 
~dl. 
du4 
~d£ 
dug 
~dl 



— = 2r 



£1*4 



3«4 

1 + r 



9wf 



(1+r) 2 
-2(l-e)ue- 



10ue 
1 + r 

45^4^6 



27m 4 3 



(1 + r) 2 (l + r) s 



(A2a) 
(A2b) 
,(A2c) 



where e = 4 — d. Linearization of these flow equa- 
tions about the Wilson-Fisher fixed point (r* ,u\,u1) = 
(-e/6 + (9(e 2 ), e/9 + 0(e 2 ), e 3 /54+0(e 4 )) yields three RG 
eigenvalues 



Ai = 2-e/3 + C>(e 2 ) 
A 2 = -e + 0(e 2 ) , 
A 3 = -2-3e + 0(e 2 ) 



(A3a) 
(A3b) 
(A3c) 



Ai = \Jv > Ois the inverse correlation length ex- 
ponent, A 2 is the scale dimension of the least irrelevant 
operator, and A3 is irrelevant even for e = 0. The in- 
frared properties of the triangle diagram thus do not 
lead to another (besides A2) eigenvalue of 0(e). This is 
true even though it does lead to a leading scaling be- 
havior of uq that is given by UQ(b — + 00) ~ b~ e , as 
can be seen by solving the flow equations. However, 
this just means that u 6 is not a proper scaling opera- 
tor, since it couples to 114. With Su^ — u± — u% and 



8u 6 = u 6 



Hi ; 



the proper next-to-least irrelevant scaling 



operator is g 6 = Su 6 — (e 2 /2)<5u 4 + 0(e 3 ). Its scale di- 
mension is [ge] = A3 , so gg is indeed irrelevant with scale 
dimension —2 in d = 4. 



This phenomcmon of the generation of a new oper- 
ator that is irrelevant with respect to the perturbative 
fixed point in d = 4 — e is completely analogous to the 
irrelevance of the three-body interaction with respect to 
the perturbative metal-insulator transition fixed points in 
d = 2+e, and the analogy extends to the long-rangedness 
of the new interaction. It is interesting to note, however, 
that there is no general a priori guarantee that ge will no 
be relevant in, say, d = 3. Indeed, the bare scale dimen- 
sion of i*6 is —2 + 2e, which naively implies a marginal 
operator in d = 3. The one-loop correction switches the 
sign of the O(e) correction, see Eq. (|A3c|l . so the one- 
loop approximation to A6 makes g$ more irrelevant with 
increasing e, not less. However, it is important to keep 
in mind that, (1) this change of sign is a special prop- 
erty of </> 4 -theory, which can be seen only by means of an 
explicit calculation, and, (2) there is no guarantee that 
high-order terms in the e-expansion will not have a pos- 
itive sign and make A3 positive in d — 3. In ^ 4 -theory 
there are no indications that this is the case, but it could 
happen in a different, more complicated field theory. 
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Even though there are qualitative differences between an 
Anderson-Mott transition and an Anderson transition, 
e.g., the former has a critical density of states, this state- 
ment is true in the following sense: Both transitions are 
driven by diffusive modes, and the lower critical dimen- 
sionality is equal to 2. The presence of the interaction just 
modifies the diffusive modes. A Mott transition, by con- 
trast, with or without disorder, is fundamentally driven by 
interactions. 



